Complete Electric Dipole Response and the Neutron Skin in ^osp^-, 
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A benchmark experiment on ■^°®Pb shows that polarized proton inelastic scattering at very forward 
angles including 0° is a powerful tool for high-resolution studies of electric dipole {El) and spin 
magnetic dipole (Ml) modes in nuclei over a broad excitation energy range to test up-to-date nuclear 
models. The extracted El polarizability leads to a neutron skin thickness rskin = 0.1561q q2i f™ 
in ^°^Pb derived within a mean-field model [Phys. Rev. C 81, 051303 (2010)], thereby constraining 
the symmetry energy and its density dependence, relevant to the description of neutron stars. 



PACS numbers: 25.40.Ep, 21.10.Re, 21.60.Jz, 27.80.- 

The electric dipole {El) response of nuclei is domi- 
nated by the giant dipole resonance (GDR), a highly ex- 
cited collective mode above the particle emission thresh- 
old [Tl. Its properties are well understood but recent in- 
terest focusses on evidence for a soft mode in neutron-rich 
nuclei below the GDR termed pygmy dipole resonance 
(PDR). Because of the saturation of nuclear density, ex- 
cess neutrons might form a skin whose oscillations against 
an isospin-saturated core should give rise to a low-energy 
El mode 3 • Therefore, the PDR may shed light onto the 
formation of neutron skins in nuclei [3[ . Another quan- 
tity related to nuclear El modes is the symmetry energy 
acting as restoring force. The El strength distribution 
carries information on its poorly known magnitude and 
density dependence indispensable ingredients for the 
modeling of the equilibrium properties of neutron stars 

i- 

A case of special interest is the doubly magic nucleus 
^"^Pb. In a measurement of parity-violating elastic elej> 
tron scattering at JLAB, the PREX collaboration |(|| 
aimed at the first model-independent determination of 



the neutron skin thickness in ^°^Pb. However, the re- 
cent result r-gkin — 0.34lg;i[y fm suffers still from limited 
statistics. Studies of energy density functionals (EDFs) 
7 1 using Skyrme forces [8| or a relativistic framework 
9| suggest the nuclear dipole polarizability ao as an al- 
ternative observable constraining both neutron skin and 
symmetry energy. The polarizability is related to the 
photoabsorption cross section cjabs [lOj l 
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where w denotes the photon energy. Because of the in- 
verse energy weighting, ajj depends on the El strength 
at low energies. Theoretically, advanced methods exist 
in closed-shell nuclei for a realistic description of the El 
strength distributions ll|. 

The centroid of the PDR lies typically in the vicinity of 
the neutron emission threshold (S'„). Data on the PDR 



in very neutron- rich nuclei are still scarce jl2l - ll4l | . Stable 
nuclei at different shell closures have been explored with 
the (7,7') reaction (Ref. [lH and refs. therein). While 
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this technique provides important information on the fine 
structure of the PDR, it is essentiaUy hmited to excita- 
tion energies up to S'„, and unobserved branching ratios 
of the 7 decay to excited states may require corrections 
for the total strength [lej]. Measurements of decay neu- 
trons are constrained to energies > Sn and uncertain- 
ties in the vicinity of Sn are large. We present here a new 
experimental tool, viz. polarized proton scattering at an- 
gles close to and including 0°, to provide the complete El 
response in nuclei up to excitation energies well above the 
region of the GDR. At proton energies of 200 - 400 MeV 
the cross sections at small momentum transfers are dom- 
inated by isovector spinflip-Afl transitions (the analog 
of the Gamow- Teller mode) and b y C oulomb excitation 
of non-spinflip El transitions 171, ll8|. A separation of 
these two contributions, necessary for an extraction of 
the El response, is achieved by two independent meth- 
ods: a multipole decomposition analysis of the angular 
distributions (MDA) and the measurement of polariza- 
tion transfer (PT) observables. 

The ^''*Pb(p, p') experiment was performed at RCNP, 
Osaka, Japan. Details of the technique can be found in 
[igj . In the present work [i^ , a proton beam of 295 MeV 
with intensities 2 — 10 nA and an average polarization 
Pq ~ 0.7 bombarded an isotopically enriched ^°*Pb foil 
with an areal density of 5.2 mg/cm^. Data were taken 



with the Grand Raiden spectrometer [2l| in an angular 
range 0° — 2.5° and for excitation energies £^a; — 5 — 23 
MeV. Sideway (S) and longitudinally (L) polarized pro- 
ton beams were used to measure the polarization transfer 
coefficients 22| Dss' and Dll', respectively. Additional 
data with unpolarized protons were taken at angles up 
to 10°. Utilizing dispersion matching techniques, a high 
energy resolution AE ~ 25 keV (full width at half max- 
imum) could be achieved. 



sitions at low excitation energies, a resonance-like struc- 
ture close to Sn = 7.37 MeV and the prominent isovector 
giant dipole resonance (IVGDR), peaked at Er^ « 13.4 
MeV with pronounced fine structure, are observed. The 
total spin transfer S can be extracted from the mea- 
sured PT observables which at 0° takes a value of one 
for spinflip (AS" = 1) and zero for non-spinflip (AS* = 0) 
transitions Figure [D^b) shows T, lor E^ — 5 - 22 

MeV. Values between and 1 result from a summation 
over partially unresolved transitions with different spin- 
flip character. The data reveal a concentration of spinflip 
strength in the energy region 7 — 9 MeV, where the spin- 
Mi resonance in ^°®Pb is located [3], while the bump 
between 10 and 16 MeV has AS* = character consistent 
with the excitation of the GDR. The AS = 1 strength 
above the GDR may result from the AS" components of 
the onsetting quasifree scattering cross section [241] . 

A multipole decomposition was performed for angu- 
lar distributions of the cross sections in the PDR and 
GDR regions. Theoretical angular distributions were cal- 
culated with the code DWBA07 j25i] using microscopic 
quasiparticle-phonon model (QPM) wave functions [26| 
and the Love-Franey effective proton-nucleus interaction 
[27|. The interference of Coulomb and nuclear contribu- 
tions to the cross sections was taken into account for El 
transitions. For a satisfactory description of the data it 
was sufficient to include, besides El and Ml, one higher 
multipole representative for all other contributions. The 
latter was chosen to be either E2 or E3 in the region of 
the PDR. In the GDR region the Ml contribution was 
zero within error bars and was replaced by a phenomeno- 
logical background describing the data at high excitation 
energies. The weight of each component was determined 
by a least-squares fit to the data. 
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FIG. 1: (Color online) (a) Spectrum of the ^°^Ph{p,p') reac- 
tion at Ep = 295 MeV with the spectrometer placed at 0°. 
(b) Total spin transfer S. 

Figure [IJa) displays a spectrum at 0° . Strong tran- 
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FIG. 2: (Color online) Decomposition of non-spinflip {AS = 
0) and spinflip (AS = 1) cross section parts based on the 
MDA and PT, respectively, in the excitation energy region 
5 — 9 MeV. The hatched areas indicate the experimental un- 
certainties. Excellent agreement between the two completely 
independent methods is observed. 
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Cross sections for AS* = and 1 from the MDA and PT 
analysis for < 9 MeV are compared in Fig. [2] Within 
uncertainties the correspondence between the two com- 
pletely independent decomposition methods is excellent. 
This puts confidence in the MDA results discussed in the 
following, which provide much better resolution because 
of the superior statistics compared to a double scattering 
measurement of PT. In the GDR region no direct com- 
parison is possible because of the unknown AS' content of 
the phenomenological background. However, both meth- 
ods agree that AS" — 1 contributions are very small. 



tion [33| experiment. Again, very satisfactory agreement 
of all three measurements is observed. 
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FIG. 3: (Color online) (a) B{E1) strengths in '^°^Ph in the 
region Ex — 4.8 — 8.2 MeV as deduced from the present work 
in comparison with (7,7') and (71,7) experiments [2^. [29l - l3ll ] . 
(b) Photoabsorption cross sections in the GDR region from 
the present work compared to (7, xn) [s^] and total photoab- 
sorption [3^ measurements. 

Next we show that reliable B(£'l) strengths can be 
extracted from the {p,p') data. While the angular de- 
pendence of El transitions is generally state-dependent 
because of the Coulomb-nuclear interference, cross sec- 
tions at very small angles (0iab < 1°) arise purely from 
Coulomb excitation. Thus the conversion from cross sec- 
tion to strength is straightforward using semiclassical 
theory The B(£'l) distribution up to 8.2 MeV is 

compared in Fig. El^a) with an average over all available 
208pb(7,7') and 207pb(n,7) data (Refs. [H, [IMll and 
refs. therein). Excellent agreement is obtained up to Sn- 
The excess strength in the {p,p') data above the neutron 
threshold can be attributed to previously unknown neu- 
tron decay widths of the excited 1~ states, which modify 
the branching ratios in the 7-decay experiments and thus 
the extracted B(£'l) values. Figure [3ljb) shows the pho- 
toabsorption cross sections in the GDR region together 
with results from a (7, xn) fs^] and a total photoabsorp- 
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FIG. 4: Experimental B(_E1) strength distribution in ^°*Pb 
in comparison to QPM and RTBA calculations described in 
the text. Note the different scales below and above 8.2 MeV. 

Figure[4lja) displays the experimental B{E1) distribu- 
tion. From the numerous computations of the El re- 
sponse in ^"^Pb we show in Fig. |H[b) recent results from 
the QPM [2^, and (c) the relativistic time-blocking ap- 
proximation (RTBA) [11] ■ The QPM calculations con- 
tain up to 3-phonon configurations for E^ < 8.2 MeV 
and 2-phonon configurations in the GDR region. Al- 
though the RTBA has recently been extended to include 
the full set of 2-phonon states [1^ , the results shown are 
based on a particle- hole(8)phonon model space 34|. In 



the low-energy region, the QPM provides a realistic de- 
scription of the fragmentation but the overall strength 
is somewhat too small, while the RTBA model space is 
not yet sufficient to reproduce the fine structure, and the 
strength is somewhat too large. The width of the GDR 
is roughly reproduced by both models. Within the QPM 
the effective isovector interaction strength is adjusted to 
the experimental GDR centroid at 13.4 MeV. The RTBA 
calculations are fully self-consistent and the GDR cen- 
troid determined by the covariant EDF parametrization 
amounts to 12.9 MeV for the NL3 parameter set used. 
Such a comparison between high-precision data and the 
3-phonon version of the QPM guides the next genera- 
tion of self-consistent extensions of the covariant EDF. 
Taking into account higher-order configurations, ground 
state correlations, and pairing vibrations should improve 
agreement with the data. 

Finally, as discussed above, an important quantity is 
the electric dipole polarizability. We find an = 18.9(13) 
fm'^/e^ for the El strength up to 20 MeV. By taking an 
average of all available data including excitation energies 



4 



b 
>. 

15 
ra 

N 

ro 
o 
Q. 
o 
o 

Q. 

b 



25 


- ^ Exp. 






208p^ 


23 


H Ref. [8] 








21 














19 










17 







>i. 




15 












0.12 0.14 


0.16 


0.18 


0.20 0.22 



Neutron Skin Thickness (fm) 

FIG. 5: (Color online) Extraction of the neutron skin in ^°*Pb 
based on the correlation between rskin and the dipole polar- 
izability an established in Ref. @]. 



up to 130 MeV [32, l33[, a result with further reduced 
uncertainty ao = 20.1(6) fm'^/e^ is obtained. The co- 
variance ellipsoid of the correlation between a d and the 
neutron skin thickness rskin in the approach of Ref. Q is 
shown in Fig. [5l Only with the present precision for 
(hatched band) one can constrain the neutron skin thick- 
ness to Tskin = 0.156^g Q2^ fm. The hitherto most precise 
determinations of this quantity for ^"^Pb [s^fs?'] deduced 
from exotic atoms (rs^in = 0.18 ± 0.02 fm) and hadron 
scattering (rgkin — 0.21ltQ Qg3 fm), respectively, are in 
excellent agreement with our result based on a totally in- 
dependent method. Recent calculations of neutron mat- 
ter and neutron star properties [39] in the framework of 
chiral effective field theory suggest rskin = 0.17±0.03 fm. 
The predictions are sensitive to three-nucleon forces, 
which may be further constrained by the present results. 
Since the correlation between polarizability, neutron skin 
thickness and symmetry energy is model-dependent, viz. 
Tskin oc • flsym [S] , a Systematic study with a variety of 
EDFs as well as experimental tests in other nuclei would 
be important. 

To summarize, polarized proton scattering at very for- 
ward angles is a tool to study, with high resolution, the 
complete electric dipole response of nuclei from low exci- 
tation energies up to the GDR. The El strength distri- 
bution deduced in a benchmark experiment on ^"^^Pb is 
in excellent agreement with available data. It provides, 
however, new information in the region around the neu- 
tron emission threshold where all previous experiments 
had limited accuracy. A precise value for the El polariz- 
ability can be extracted with important consequences for 
a determination of the neutron skin and the symmetry 
energy in neutron-rich nuclei. Although controversially 
discussed r-skin may independently be derived from a 
similar correlation with the PDR strength [l^l , which 
is accurately determined by the present data as well. Be- 
yond these results, the experiment also confirms the spin- 
Mi resonance in ^"^^Pb. Furthermore, the fine structure 
of the dipole modes contains information on level densi- 



ties [40[ and characteristic scales [4l|, giving insight into 
their dominant damping mechanisms. 
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